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Abstract. Light degrees of freedom that modify gravity on cosmological scales must be 
“screened” on solar system scales in order to be compatible with data. The Vainshtein 
mechanism achieves this through a breakdown of classical perturbation theory, as large in¬ 
teractions involving new degrees of freedom become important below the so-called Vainshtein 
radius. We begin to develop an extension of the Parameterized Post-Newtonian (PPN) for¬ 
malism that is able to handle Vainshteinian corrections. We argue that theories with a 
unique Vainshtein scale must be expanded using two small parameters. In this Parame¬ 
terized Post-Newtonian-Vainshteinian (PPNV) expansion, the primary expansion parameter 
that controls the PPN order is, as usual, the velocity v. The secondary expansion parameter, 
a, controls the strength of the Vainshteinian correction and is a theory-specific combination 
of the Schwarzschild radius and the Vainshtein radius of the source that is independent of 
its mass. We present the general framework and apply it to Cubic Galileon theory both 
inside and outside the Vainshtein radius. The PPNV framework can be used to determine 
the compatibility of such theories with solar system and other strong-field data. 
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1 Introduction 

General Relativity (GR) is a very successful theory. Since its inception, it has not only led to 
the prediction of new physical phenomena, but has also successfully passed all experimental 
tests to this day [1, 2], But the question of whether GR is the correct theory for describing 
all gravitational phenomena in Nature still stands. Indeed, experimental gravitation is a very 
active field in physics, aimed at increasing the precision of gravitational tests, but also at 
pushing the boundaries where these tests have been carried out. 

In order to test gravity, one may simply take a gravitational theory, calculate its pre¬ 
dictions for the tests of interest, and compare with observations. However, as the “space” 
of gravitational theories is vast, perhaps infinite, this is not a very economical process. A 
better way is to construct a “master theory”, a framework, that encompasses a wide range 
of theories. Different theories would then correspond to a specific set of parameters of this 
framework. Examples are the Parametrized Post-Newtonian (PPN) [3-5], the Parametrized 
Post-Keplerian (PPK) [6], and the Parametrized Post-Einsteinian (PPE) [7-9] formalisms. In 
the cosmological regime we also have the various simplified parametrized approaches [10, 11], 
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the small-scale parametrized approach [12] as well as more complete frameworks such as 
the Hu-Sawicki Parametrized Post-Friedmannian formalism [13, 14], the Parametrized Post- 
Friedmannian (PPF) formalism [15-17], the effective fluid [18, 19] and effective field theory 
formalisms [20-22], 1 

Tests of gravity in the solar system have reached incredible precision. With the help 
of the PPN formalism, solar system tests such as lunar laser ranging and doppler tracking 
of the Cassini spacecraft put bounds on the PPN parameters of around 10~ 3 — 10 - ' 5 for 
curvature effects and 10 -7 — 10~ 20 for preferred frame effects [2], However, the PPN for¬ 
malism (as well as the PPE formalism) has as a basic tenet that the spacetime away from a 
source is asymptotically flat and the gravitational theory remains perturbative down to the 
Schwarzschild radius of the source. But there are theories for which the latter assumption is 
manifestly untrue. Those theories have the property that the additional degrees of freedom 
to the metric become strongly coupled at some macroscopic scale, which in turn makes it 
impossible to construct a perturbative expansion (as the PPN requires) that would be valid 
from infinity all the way down to the Schwarzschild radius. But why should we bother with 
such theories? 

The theories in question were constructed as a way to introduce departures from GR in 
the cosmological regime, i.e. in the limit of ultra-low curvatures and potentials. Whilst it is 
of academic interest to explore this possibility on purely theoretical grounds, most of these 
theories seek to address the cosmological constant problem [23] and/or the origin of cosmic 
acceleration through a modification of gravity (for a review, see [24]). Cosmic acceleration 
has by now been substantiated by a wide variety of observations. Although hints for cosmic 
acceleration may be traced in studies of large scale structure in the 80’s, the first real evidence 
came from measurements of the luminosity distance of type la supernovae [25, 26]. Cosmic 
acceleration is thought to be caused by an effective fluid, called Dark Energy, whose nature is 
still a mystery. The latest supernovae results from a combination of Sloan Digital Sky Survey 
(SDSS-II) and Supernovae Legacy Survey (SNLS) data constrain the equation of state of 
Dark Energy to w = P/p ~ —1.018 ± 0.057 [27], where P is the dark energy pressure and p 
its energy density. Other observational probes also indicate a Dark Energy component, for 
instance, the latest measurements of the Cosmic Microwave Background (CMB) anisotropies 
from the Planck Surveyor [28] give w = —1.13)]q']q while the cross-correlation between the 
Integrated Sachs-Wolfe (ISW) effect from the CMB and Large Scale Structure (LSS) favour 
w = — 1.0ltg4Q, at around 4<r [29]. 

Attempts to explain cosmic acceleration by modifying the theory of gravity on cosmo¬ 
logical scales ultimately require new degrees of freedom that make gravity behave rather 
differently from GR on ultra low curvatures. However, since the solar system data indicate 
that gravity is described by GR to a very good approximation in this regime, these degrees of 
freedom must somehow be hidden there. As solar system data (but also other data towards 
the strong gravitational field regime) are probing curvature scales much larger than cosmolog¬ 
ical curvatures this is entirely possible, at least observationally [30]. Indeed, this can happen 
if classical perturbation theory involving the new degrees of freedom breaks down at some 
large distance from a massive source, far beyond the Schwarzschild radius. In the case of the 
Sun, we typically require this large distance scale, known as the Vainshtein radius, to extend 
beyond the edge of the solar system, so that “screening” occurs within the solar system itself 
thanks to the nonlinear physics. The mechanism by which this occurs is called the Vainshtein 

1 Although [13] and [17] use the same acronym, the two formalisms are widely different. 
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mechanism [31-33], and is exploited by so-called Galileon theories [34-36], massive gravity 
[37, 38], k-Mouflage gravity [39], and the Fab Four [40-42], to name just a few. An effective 
field theory for the Vainshtein mechanism has also been developed [43] using the Horndeski 
action [44, 45]. However, the screening of these new degrees of freedom is not perfect, and 
some residuals can in principle be detected. This paper presents an extension of the PPN 
framework, the Parametrised Post-Newtonian Vainshteinian (PPNV) framework that is able 
to handle Vainshteinian corrections and thus paving the way for determining the compati¬ 
bility of such theories with solar system and other strong-field data. There has been some 
previous work in this direction, most notably [46-50] . Although it is easy to imagine how one 
can build a PPNV framework outside of the Vainshtein radius where classical perturbation 
theory remains valid, it is less obvious to see how to do this in the interior where the clas¬ 
sical perturbative description has broken down. Salvation lies in the so-called classical dual 
description [51, 52] , which ultimately corresponds to a Legendre transform of the original 
theory and admits a classical perturbative expansion inside the Vainshtein radiius (but not 
outside). We will exploit this procedure in developing our framework. This initial work will 
focus on theories for which there is a unique macroscopic scale beyond the Schwarzschild 
radius at which point classical non-linearities start to kick in. 

Whilst the need to develop a PPNV formalism is clearly important in testing these 
particular theories, we ought to issue a gentle word of warning about their validity once the 
non-lineariities begin to kick in. The point is that the break down of classical perturbation 
theory is inherited from a breakdown in perturbative unitarity above some mass scale, A. 
From an effective field theory perspective, one would certainly expect there to be a tower of 
additional higher dimensional operators suppressed by the same scale A, included to preserve 
unitarity. As emphasised in [53], such operators could, in principle, affect the classical solution 
out to macroscopic scales, possibly up to the Vainshtein radius itself. Without knowing the 
details of the (partial) UV completion beyond A, we do not know the impact of these operators 
but in principle they could affect the details of the PPNV expansions for a given theory 
both outside and inside the Vainshtein radius. We emphasise that this paper seeks to build 
towards a PPNV framework for generic theories, which ultimately may be taken to include 
the aforementioned higher dimensional operators. However, for simplicity of presentation 
we focus on the Cubic Galileon [34] as an illustrative example, and do not worry about the 
precise details of which operators are present in that case, and that might have been omitted 
from a consistent UV completion. 

The paper is organised as follows. We first describe the PPN framework, followed by 
a short description of the Vainshtein mechanism. We then outline how the PPN formalism 
has to be modified in order to include Vainshteinian corrections. As we will show further 
below, there are two separate perturbative regimes, one inside the Vainshtein radius and one 
outside, which must be handled separately. We then take the Cubic Galileon theory as an 
example and perform the calculation in the two regimes. We finish by considering spherical 
symmetry as a special case. 

2 Vainshteinian corrections to the PPN formalism: the PPNV expansion 
2.1 The Standard Post-Newtonian Approximation 

The PPN formalism is a prescription for a perturbative expansion of the gravitational, matter 
and additional field equations of motion in successive orders of a small parameter, the velocity 
of matter v (in units of the speed of light). It was developed by Nordtvedt[3] and later 
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expanded refined by Will [4] and Will and Nordtvedt [5] . Here we review the basic ingredients 
of PPN. 

In PPN, one assumes that gravity is a geometric phenomenon and is described by a 
metric g^ v and that matter only “sees” this metric g^ and follows its geodesics. Specifically 
this means that if T^ v is the stress-energy tensor for matter then 

v /t n = o, ( 2 .i) 

where is the covariant derivative compatible with g^ u . However, the metric g^ may not 
be the only gravitational field. For instance, there may be additional scalar fields to the 
metric that mediate the gravitational force (as is in our case), however, these fields do not 
enter the equations of motion for matter when those are formulated in the Jordan frame. 
Indeed, it is the so-called Jordan frame metric that should always be used in PPN. 

To proceed, one first expands the metric g liv as a weak field perturbation h^ u around 
Minkowski spacetime as 

dflV = T hjxw (2-2) 

When doing explicit calculations, we shall adopt the convention that the metric has signature 
(—, +, +, +). The metric perturbation h /iu is further expanded in successive orders dictated 
by a small parameter: the velocity of matter v (we will assume units where the speed of light 
is unity). Thus, the PPN is a small velocity expansion and is valid for gravitational fields 
generated by slowly moving matter. 

The first correction to Minkowski spacetime is the Newtonian potential 

U(x) = [ P( " X ’ t ] i (fix' ^ V * 2 U = —4t rp, (2.3) 

J |x-x'| 

generated by the rest-mass density of matter p. This potential generates acceleration a = 
GatVH where Gat is the measured Newton’s constant by a lab experiment on the Earth. 
Thus, virial relations dictate that v 2 < G^U, implying that the Newtonian potential has a 
PPN order 0(2). In particular for spherical symmetry, we have that 2GnU = (f < 1 where 
r s = 2 GnM is the Schwarzschild radius of a source with mass M. Moreover the acceleration 
equation gives the PPN order for derivatives, i.e. spatial derivatives leave the PPN order 
unchanged while time derivatives have a PPN order ^ ~ u • V ~ CJ( 1). Of course, these 
considerations are not strict equalities and should be taken as a rule-of-thumb in order to 
assign PPN orders to the various quantities, such as v or U that may enter the equations of 
motion. 2 

Other matter variables are also associated with a PPN order. The matter rest-mass 
density p has a PPN order 0(2) as is implied by (2.3). In addition to the density, matter 
sources may also have specific energy density n (coming from other forms of energy, e.g. 
compressional, thermal, and corresponding to the ratio of the total of these energy densities 
and the rest mass density of matter p) and pressure P, so that the matter stress-energy 
tensor takes the form 

Tf lI/ = (p + pU + P)Uf i u l/ + Pg tlv , (2.4) 

2 Spherically symmetric orbits around a point mass imply strict equality v 2 = GnU while a hypothetical 

2 

body oscillating through the center of the Sun would have v RM $ ~ -nP — GnUo where r o is the extend of 

h sun 

its motion and Uo is the gravitational potential at the centre. Hence once again v RMS < GnU as ro can be 
chosen at will. 
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where the velocity four-vector is normalized as g flu ,u fJ ’u u = — 1. Since the pressure is generally 
smaller than the gravitational energy of matter ( P < pU ) we associate it with a PPN order 
P ~ 0(4). Similar considerations apply to the specific energy density II < G^U that is 
associated with a PPN order II ~ 0(2). The evolution of the matter fluid is determined by 
the continuity and Euler equations of motion that are obtained from (2.1) by expanding in 
PPN orders. 

Clearly then, (2.4) implies that 0(4) terms would appear in the (generalized) Einstein 
equations, thus we need to include these terms in the metric perturbation. The required order 
for the metric perturbation components is found to be hoo = where the superscript 

number in the parenthesis denotes the PPN order of a particular term, hoi ~ 0(3) and 
hij ~ 0(2) [54]. In the standard PPN formalism one has that /?q 0 = 2 G^U, however, as we 
shall see further below, in the PPNV framework the Vainshteinian corrections contribute to 
/ig 0 in the form of new Vainshteinian potentials. Furthermore, in the standard PPN gauge, 
one has h t] = ‘2G njU 5jj where 7 is one of the PPN parameters. In the PPNV example we 
consider further below, one still has that h tJ oc 5ij, however, h t] is no longer proportional to 
hgg' as in PPN due to the appearance of new Vainshteinian potentials. We shall show that 
PPNV introduces new parameters that are not part of the standard PPN parameters. 3 

In general, in PPN, one proceeds by determining the Too equation up-to 0(4), the 
Roi equation up-to 0(3) and the Rij equation up-to 0(2). Solving these equations in turn 
determines the metric to the required order and provides the definition of the PPN potentials. 
In addition to U, in the standard PPN gauge, there are a further 8 PPN potentials that appear 
in orders higher than 2 . Since in this article we are mostly concerned with the presentation 
of the PPNV extension compared to the PPN formalism, we shall not consider the 0(3) and 
0(4) PPN orders, which are left to future work. Therefore we will only need the I?oo and 
Rij equations and the stress-energy tensor to 0(2). These are 

B^ = -\v 2 hg\ (2.5) 


and 


Rij — 0 V A 


Vih k j + Vjh k i 




( 2 . 6 ) 


where h = h l t . We will also need Too = P and the trace T = T^ = —p. Let us now 
discuss the Vainshtein mechanism and the modifications it introduces, leading to the PPNV 
expansion. 


2.2 Vainshtein Mechanism for a Spherically Symmetric Source 

As we already discussed in the previous subsection, the metric theories treated with the stan¬ 
dard PPN formalism (e.g. GR or Brans-Dicke [54]) have a single expansion parameter given 
by v 2 In contrast, the non-linearities of the scalar field in theories with a Vainshtein 

mechanism (such as the Galileon theories) introduce new regimes that must be handled along 
with PPN. These Vainshteinian thresholds are associated with an extra expansion parameter 
^7 that accounts for new physical effects associated with the new scale introduced by these 
theories, the Vainshtein scale ry. 

A simple way to think about their different expansions is the following. The ratio 
y < 1 accounts for the non-relativistic nature of the matter sources considered in the PPN 

3 This is not a unique feature of the PPNV formalism. For instance in general bi-gravity theories a similar 
situation arises [55]. 
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approximation while the ratio determines whether the observer is placed either inside or 
outside the Vainshtein radius ry and how far from the source (or close to ry ) and consequently 
how important the Vainshteinian corrections are. Clearly, we must then consider two regimes: 
inside the Vainshtein radius so that — <C 1 and outside the Vainshtein radius where — S> 1. 

Ty Ty 

These two regimes must be treated separately. 

As an example consider static spherically symmetric solutions in the Cubic Galileon 
theory in the Einstein frame (see Eq. 3.1). The metric that enters (2.1) is a combination 
of the Einstein metric g flu and the scalar field x (assumed dimensionless) as g^ v = e^g^. 
Thus, considering weak fields on Minkowski spacetime the effective gravitational potential 
would be /ioo = ^oo — 2y = 2G f M — 2%. The scalar equation of motion for a spherically 
symmetric source of mass M centered at r = 0 is 


2cj + 3 d r r 2 n , r r (v^2l = 

r 2 I r ^ J + r 2 rfry ) \ ^2 ’ 


r 2 dr 


r* 


where 


a = 


Mp 
A 3 ' 


Now (2.7) can be integrated once to give 

(2a; + 3)r 2 y' + ar(x') 2 = GM. 

Thus, solving for x! we find 


. 2a;+ 3 

v = - r 

2 a 


—1 + \ 1 + 


4 GMa 


(2a; + 3) 2 r 3 


(2.7) 

( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


where the plus sign in front of the square root is chosen so that no divergence appears when 
a —> 0. There are two limiting cases. As a —>■ 0, or in other words, when ( 9 ^+ 3)? r 3 ^ 1 then 

GM (GM) 2 a 

X “ ~~ (2w + 3)r + 4(2 cj + 3 ) 3 r* + " ' ^ 


In the opposite limit asa-l 00 , or in other words, when ( 2 t+ 3 )V 3 ^ 1 then 


X = 



( 2 . 12 ) 


where the integration constant is ignored as it can be rescaled away by a coordinate trans¬ 
formation. 

The turning point occurs when = 1, he. at a radius ( 2 7 r) 1 / 3 ( 3 + 2 a ;) 2 / 3 ry w ^ iere 

ry is the scale given by 


ry = 


A 


M 


M, 


pi 


1/3 


(2.13) 


This ry is called the Vainshtein scale and will play a fundamental role in what follows. 

The Vainshtein mechanism is now clear. For small distances away from the source, 
r<ry, the effective gravitational potential is 


hoo 


r s 


r 




(2.14) 
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while for large distances away from the source, r ry, it is 

1 /ry \ 3 


hoo = ~ 
r 


1 - 


64tt(2o; +3) 2 (2+ w) \ r J 


+ 


(2.15) 


where r s = j^r,r s is a rescaled Schwarschild radius. Note that with ry now defined, the a 

r 3 

parameter can be written as a ~ y- 

We saw by example what we should expect in a theory that leads to a Vainshtein 
mechanism. In the next subsection we formalise the treatment above and consider general 
theories with a Vainshtein mechanism. In particular we shall see how the Vainshtein scale and 
the expansion parameter a emerge, and what impact this can have on the PPNV expansion. 

2.3 The PPNV formalism 

To develop a Parametrised Post-Newtonian Vainshteinian (PPNV) formalism, we must first 
understand some of the generic features of modified gravity theories with Vainshtein screen¬ 
ing. The modification of gravity is normally associated with a light scalar degree of freedom 
that strongly affects the dynamics on cosmological scales when linearised perturbation theory 
is valid. However, in the presence of a heavy source, linearised perturbation theory breaks 
down at some macroscopic scale owing to derivative interactions involving the new field y. 
Let us illustrate this with a generic example, schematically described by an action in the 
Einstein frame 4 : 


S 


d A x 


Mpy/—gR - (. d x f + 0(x) + Mph^O^ix) + + tt? 


M rt 


(2.16) 


where M 2 = and where in general G/Gn ~ 0(1) (usually the measured Newtonian 
strength Gat is not equal to the bare gravitational strength G in the action but is typically 
related to it by an 0(1) quantity). This schematic form encompasses a number of modified 
gravity scenarios in the so-called decoupling limit, including single field Galileons [34] and 
dRGT massive gravity 5 [37, 38]. The canonical scalar field is coupled directly to matter 
with gravitational strength, but the Vainshtein mechanism is made possible by the derivative 
interactions between the scalar and the graviton and/or the self interactions. Assuming for 
the moment a unique strong coupling scale, A, we have that these operators contain terms 
like 


0(x) 3 


d 2r 


'X 


0^{x) a 


d 2n 


'X 


n— 1 


(2.17) 


^2m-\-n —4 ’ J^2m-\-n —4 ’ 

where the integers m, fh > 1, n, n > 3. The strong coupling scale marks the breakdown 
of perturbative unitarity, and in principle one may have to include a whole tower of higher 
dimensional operators to preserve unitarity at higher energies. 

For a static and spherically symmetric profile, the classical potentials at large distances 
from a heavy source of mass M are simply Newtonian, M p h c ^ u ~ Xu n ~ • Assuming this 

Newtonian behaviour, the interactions in (2.17) would become comparable with the leading 
order canonical terms when 


1 


1 


j^2m-\-n—4 y>2m 


M 1 

M p r 


1 ( M 1 

M p r 


J* r-u - 


M 

w n 


2m-\-n—4 


4 If gravity is modified by a long range scalar, you would expect the leading order canonical kinetic structure 
presented here in any sort of sensible Lorentz invariant set up that is weakly coupled at long distances. 

5 In massive gravity, the scalar y is identified with the Stuckelberg scalar.The vector modes are subleading 
in the decoupling limit, so our method would have to be adapted to include them. 
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with similar expressions for the barred integers. We see that in general different interactions 
yield different macroscopic scales. As a first pass in developing the post-Vainshteinian for¬ 
malism, let us only consider scenarios in which there is a unique macroscopic scale beyond 
the Schwarzschild radius, so that all nonlinear interactions become important at the same 
scale. In other words, we only consider those theories for which there is a unique strong 
coupling scale, A, and for which all interactions of the form (2.17) yield the same value of 
2m+n-4 = S ’ 2rh+n-i = s • An exam ple of this would be a combination of Galileon interac¬ 
tions suppressed by the same strong coupling scale, each giving s = Therefore, we shall 
begin to develop our formalism only for those theories for which there is a unique macroscopic 
scale of interest, which we then identify with a unique Vainshtein radius, 


r v ~ 


1 

A 



n — 2 

2m + n — 4 


n — 2 


2 m + n — 4 


;Vra ,m G 0(y),Vn ,m G O flv (x), 


(2 - 18) 

signaling the breakdown of classical perturbation theory. This is how the strong coupling 
manifests itself classically. Generalisations that take care of multiple macroscopic scales ( 
i.e. interactions with different strong coupling scales, and different values of s ) and multiple 
expansion parameters will be left for future study. 

Let us recall that in the PPN formalism the leading order contribution is the Newtonian 
potential, satisfying V 2 C/ = —AitGnp and is defined to be 0(2) in PPN. The orders of 
velocities and time derivatives are then inferred using the virial (v 2 ~ U) and Euler (dt ~ v- V) 
relations. In the spherically symmetric scenario described above, we note that L/ ~ —, where 
r s ~ M/M 2 is the Schwarzschild radius of the source, so there is a sense in which the PPN 
expansion is an expansion in Similarly, we shall think of PPV as an expansion in 

3 

^^7 j Although our analysis here is generic, the exponent | is ultimately motivated by 
the Cubic Galileon scenario. This is a sensible starting point since the Cubic Galileon is the 
lowest dimensional operator corresponding to a purely derivative self-interaction of a single 
scalar in a Lorentz invariant theory. 

To gain further insight, consider what happens deep inside the Vainshtein radius. If the 
scalar dynamics is dominated by a pure interaction term, we would have 


f) 2m ( v c 1 
VAnon/m/ 

^2m+n-4 


n—1 


T 

M ~ P ’ 


for some m, n, and so for a static spherically symmetric configuration, 


Xnonlin „ (1_\ (Zi 

Mp \ry) V r 


This suggests that the scalar part of the physical metric that couples to matter is 0(2) in PPN 
and O ^| 2m n--T 4 )i n PPV. However, if the deep Vainshteinian behaviour is dominated by 
a mixed interaction then we would schematically have 


Mphn V 


Q 2 ™(X 


c \n—2 

nonlin) 


^2m+n—A 


T 

Mp' 
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for some m, n. For a static spherically symmetric configuration in which screening is effective, 
we expect M p d 2 h llv ~ and so 


Xnonlin ( r \ ^ ( r _s 

M p \ry) \r 


This now suggests that the scalar part of the physical metric that couples to matter is 0(2) 


in PPN and O ( § 


2m+n—4 
n —2 


in PPV. 


The previous paragraph indicates two important things: (i) that even in the deep Vain- 
shteinian region the leading order behaviour of the physical metric is Newtonian, as it should 
be; and (ii) that what we really have here is a double expansion, owing to the hierarchical 
difference between the the Vainshtein radius and the Schwarzschild radius. To this end we 
assign a combined order 0(N, V) to a quantity that is order N in PPN, and order V in PPV. 
In terms of Schwarzschild and Vainshtein radii, we may think of 


0(N,V) ~ - 
r 


r \ 

rv) 


3V 

2 


Thus, U and xf in are 0(2,0), whilst Xnonlin is ° 2 


V 3 


2m+n —4 
n—1 


or O 


( 2 - 


2m+n—4 
n— 2 


depending on whether it is a pure scalar, or a mixed interaction that dominates the deep 
Vainshteinian region. 

PPN can be thought of as an expansion in the Schwarzschild radius, and is therefore 
equivalent to an expansion in the source. Similarly PPNV can be thought of as an expansion 
in Schwarzschild and Vainshtein radii, and for practical purposes this is best realised in terms 
of an expansion in the source, and an operator, a, that is independent of the source, but 
which carries PPV order 6 . Given the Vainshtein radius (2.18), there is only one candidate 
for a (or powers thereof), namely: 


a ~ M p / A 1 /'. 


(2.19) 


It is easy to see that a is O (—2,—^). Operators containing powers of A should now be 
rewritten in terms of a. We will also assume that velocities are always determined by the virial 
relation with the Newtonian potential, since the scalar potential never dominates over the 
(Einstein frame) graviton if screening is active. The rules of the game are therefore as follows: 
the orders for velocities (v ~ 0(1,0)), time derivatives (dt ~ 0(1,0)), space derivatives 
di ~ 0(0,0), total energy density p ~ 0(2,0), specific energy density (II ~ 0(2,0)), and 
pressure P ~ 0(4,0) are inherited from PPN [54], and do not care about the PPV expansion. 
The only operator to carry PPV order now is a ~ 0 (—2, — ^). 

For a given theory, we write the action explicitly in terms of a, and study the resulting 
field equations order by order in the double expansion. Schematically, for our generic example, 
we have 


5 ~ y d A x {M 2 [V^gR - (dxf + P(x) + V^(x)] + K- T ' 1V + XT} , (2.20) 

6 Tliis perspective is a bit sloppy in both PPN and PPNV: the source, p, is not dimensionless and therefore 
one might consider it a poor expansion parameter. In PPNV, the same may be said of a given by (2.19). 
However, we remind the reader that what is really going on is an expansion in r s /r and r/rv , with p and a 
simply being a more convenient way of keeping track of that expansion. 
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where we have rescaled x M p x, and 


P(X) = 


Q(M pX ) 

Ml 


D a n ~ 2 d 2m X n , 


P^ix) = 


0^(M p x) 

Mr, 


D a n ~ 2 d 2fh X n ~ 1 . 


( 2 . 21 ) 


Recall that we have assumed each interaction has the unique value of s = 2 m+n-A > 2 m+n-A 
where n, m E 0( X ),n,fh E 0^ v {x). The PPNV analysis is straightforward beyond the 
Vainshtein radius where the leading order contributions are purely Newtonian i.e. zeroth 
order in PPV, and corrections carry negative PPV order as higher order terms in a begin to 
kick in. These corrections generate new classes of PPV potentials, and a complete formalism 
should ultimately include all possibilities, or at the very least, all possibilities that are relevant 
to models that appear in the literature. 

In contrast, deep inside the Vainshtein radius, the leading behaviour should again be 
Newtonian but with corrections carrying positive PPV order, requiring terms proportional to 
inverse powers of a. However, the action expressed in (2.20) is not well suited to an expansion 
in terms of these inverse powers. This is where the so-called classical dual [51, 52] comes 
in. For our schematic action (2.20), there exists a dual action describing the same classical 
physics [52] 


S- 


/ ( T max 

Ml } ^gR - (d X f + £ 


r=0 L 


dF 


dA 


+ h a p 


Ml • "lSr 


dF a P 
dA 


Ml , 




-\-Qi 


—st 


F- 


E 


dF 


-A. 


A " "Ml •••Mr* 


+ O' st h a p 


^ _ g jpm a 


sa 


Ml 


Ml 


+K V T^+ X T, 


( 2 . 22 ) 


where r max is the largest value of 2m, and t = max j 2m ^ 1 4 , 2Tn M% 4 ; n, m E 0( X ),n, m E O mi/ (x)|- 
We also have combinations of polynomials in the auxiliary fields A pi __n r , i.e. 

FDF n (A,A p ,...), r^F^fi,!,,...), (2.23) 

where F n is order n in the auxiliary fields, while FS^ is order (h — 1). These polynomials 
stem from the interactions given in (2.21). Introducing 

S(m, n) = 2 [ J* + - (t - 2m n + _" 1 > 0 Vn,meO(x), (2.24) 

fi) = 2 m"+B - 4 (* - 2m _ + _" 2 ~ 4 ) >0 Vfi,m € 0<"(x), (2.25) 

it can be shown that F n scales as a~ s ^ n ~ 2 \ while F^l { scales as q,— 2 )(e-t-i) \y e now 
have an action made up of non-negative powers of 1/a, and as such it is well suited to an 
expansion that increases both the PPV order (as desired) and the PPN order. The PPV 
expansion deep inside the Vainshtein radius generates yet another class of PPV potentials, 
and due to the non-linear nature of the problem to leading order, we do not expect to be 
able to write these in closed form. 

Let us now summarize the steps that need to be followed in order to apply the PPNV 
expansion for the class of theories under consideration: i.e. those for which all interactions 
kick in at the same macroscopic scale. 
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• Determine the Vainshtein radius ry for the theory. This is usually determined from 
the Einstein frame action using (2.18) and the discussion that precedes it. 

• Form the expansion parameter a and determine its PPNV order. 

• For the outside region, reduce the field equations of the theory in the Jordan frame, as 
is done in the case of PPN. Note that a usually carries a PPN order, so care must be 
taken to take this into consideration when performing the expansion. 

• For each PPN order, solve the field equations to determine which potentials arise. 
Furthermore, for each PPN order, solve the scalar equation in PV orders to determine 
the metric potentials to the required Vainshtein order. 

• Dualize the theory in the Jordan frame and determine the held equations. 

• For the inside region, reduce the held equations of the theory in the Jordan frame, as 
is done in the case of PPN. 

• For each PPN order, solve the held equations to determine which potentials arise. 
Furthermore, for each PPN order, solve the scalar equation in PV orders to determine 
the metric potentials to the required Vainshtein order. 

We shall now develop the PPNV scheme for Cubic Galileon theory in order to illustrate the 
idea. 

3 Case study: Cubic Galileon theory 

3.1 A short introduction to Cubic Galileon theory 

In order to demostrate our formalism, we shall apply it to a particular case of a screening 
theory: Cubic Galileons. We begin by presenting the action and held equations in the 
standard form (we shall give an alternative, dual formulation further below). The action in 
the Einstein frame is 

S[g,x] = YffffQ j d ‘ ix \/-^R + j d A x^-ij (cqX + pifixj + S M [g\, (3.1) 

where G = is the bare gravitational strength, co > 0 is a constant, X = — „x 

and A is the strong coupling scale when co ~ 0(1). G and M p are not the measured value of 
Newton’s constant and Planck mass but are rather parameters in the action; we shall return 
to this point further below. 

Following the discussion above (2.18), we determine the Vainshtein scale. We have that 
for the Cubic Galileon the 0^ u term is absent while the 0{x) term is Hence, m = 2 

1 /3 

and n = 3 which gives the Vainshtein scale as ry = (yr) as i n (2.13). Having deter¬ 
mined the Vainshtein scale, we also form the expansion parameter a = The expansion 
orders for a are then found to be a ~ 0(—2, —2), in other words, a lowers both the PPN 
and the Vainshtein order of any terms multiplying it, by 2. 

The Jordan frame is determined via a conformal transformation to a new metric g pi/ = 
e 2 */ M Pg^ u that is minimally coupled to matter. Defining cf = e~ 2x ^ Mp the action in the 
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Jordan frame takes the form 


S = 


1 


16irG 


d 4 xy/^g 


2u 

<f>R + —Y 


a I j 4 ,— Y. 


- j J d 4 x^g^O(j) } + S m \g,^ A ], (3.2) 


where ui = and Y = — ^g^V n<f>V v (j). The matter action S m depends, in addition to a 
generic set of matter fields if> A , only on the metric g fW but not on the Galileon held (j). 

Varying the action with respect to g^ v and <j) and after some algebraic manipulations 
gives the generalized Einstein equations as 


fiRuu = 87 xG 


Tfiu ^ Tg^ 


uj 1 

+ v 4* + 2 


a J 

+ 8^\“ 

and scalar equation as 
'5 


YU^guv + 


Y 

6— — □</> 


Vp(/>V v (f> - VpYV v <f> — Vp(j>V„ y\, (3.3) 


a 


(3 + 2 u)0(j)+ 


V^V^V + 


18V 2 


- (D </>) 2 - -DY-Y^n<f> 


= 87 tGT. (3.4) 


We have chosen this form of the (generalized) Einstein equations as it is more convenient 
when applying the PPN (and by extension the PPNV) formalism. 


3.2 The PPNV expansion outside the Vainshtein radius 

Outside the Vainshtein radius the self-interactions of the Galileon are subdominant and 
the matter terms lead the sources of the held equations. It is straightforward to see that, 
at leading order the held equation (3.4) basically reduces to that of a Brans-Dicke held. 
In contrast to what happens inside the Vainshtein region, for sufficiently large scales the 
Galileon shows up as a linear held with small corrections due to the scalar self-interactions 
that may be responsible for modifications to the cosmological dynamics. 

Because the source is assumed to be non-relativistic it is reasonable to split the held in 
PPN orders 


^ = (l + ^ (2) + ¥> (4) + •■•), (3.5) 

where the superscript denotes the PN order of each term and <f>^ >ut ' > is the constant background 
(cosmological) value of the held outside the Vainshtein radius. We now proceed to the 
generalized Einstein equations (3.3) order-by-order. 


3.2.1 /i 00 to 0(2) 

Using (2.5) the (generalized) 00 Einstein equations outside the Vainshtein radius give at order 

0 ( 2 ) 


where we have dehned 


1 

2 


V 2 /i, 


( 2 ) 

00 


47 xGp- JvV 2) , 


G 


G 

Aout ) ‘ 

v'o 


(3.6) 

(3.7) 


12 - 














Notice that the above equation justifies ip ~ 0(2) + .... The solution to (3.6) is 

h$ = 2GU + p ( - 2 \ (3.8) 

where U is the standard Newtonian potential (2.3). 

At this point, the solution (3.8) is formally identical to the solution for /loo hi Brans- 
Dicke theory. However, we do not yet know the solution for the new potential p>^ 2 \ which 
will come out to be different than the Brans-Dicke case. 


3.2.2 hij to 0(2) 

Now consider the Rij equation to 0(2). Using (2.6) in (3.3) we find 


:V fc 




-^V% + iv i V i (/i£ ) -/i) = 4irG f r fij + ±V 2 <pW'Yi j +ViV j <pW. (3.9) 


We impose the gauge-fixing condition 


which brings (3.9) into the form 

- = ^nGfr/ij + ^V 2 ¥? (2) 7ij. 

The solution to the above equation is 


— 


2 GU - <p (2) 


7 ij- 


(3.10) 


(3.11) 


(3.12) 


With (3.8) and (3.12) we have now determined the metric to 0(2) in PPN in terms of 
the Newtonian potential U and the still unknown field p( 2 \ To proceed further we must now 
consider the scalar equation to 0(2). 

3.2.3 Scalar equation to 0(2) 

The scalar equation to 0(2) gives 

vV 2) --—^-^{(vV 2 )) 2 +v^ 2) -wV 2) -|v 2 |v^ 2 )| 2 } = --^p. (3.13) 

4(3 + 2w)4 } l 2 J 3 + 2w 

Taking a —> 0 recovers the Brans-Dicke equation for the scalar that can be easily solved. In 
the general a / 0 case, however, it is impossible to solve the above equation analytically, 
except in idealized situations, for instance, spherical symmetry. In order to make progress we 
appeal to the PPNV expansion: we can find the solution by expanding in Vainshtein orders 
as 

OO 

<p( 2 ) = <p( 2 ’°) + <p( 2 ’- 2 ) + <p( 2 ’- 4 ) + ... = J2 V (2 '~ 2n) - (3.14) 

71=0 
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3.2.4 Metric solution to 0(2,0) 


Consider first the leading Vainshtein order, i.e. 0(2,0). In this case we 
Dicke equation that is solved as 

recover the Brans- 

^ 2 ’°)= 2G U. 
y 3 + 2w 

(3.15) 

Therefore to leading order we find using (3.8) and (3.12) 


4 o 0) = 2G cU , 

(3.16) 

hf/ 0) = 2GciU lih 

(3.17) 

where we have defined the ’’cosmological” gravitational strength Gc as 


r _4 + 2cu ~ 

3 + 2cj 

(3.18) 


and recovered the PPN parameter 7 for Brans-Dicke theory: 

1 + u 


7 = 


2 + u 


(3.19) 


Thus to leading order, the Cubic Galileon theory reduces to Brans-Dicke theory outside the 
Vainshtein radius. In fact, it can be shown that this feature is valid in the PPN expanstion to 
0(4), however, due to the more complicated equations in that case we leave it for a separate 
publication. 

3.2.5 Metric solution to 0(2,—2) 

We now solve the scalar equation to 0(2, —2). Keeping on the the 0(2, —2) terms, the scalar 
equation (3.13) gives 


vV 2 ’- 2) = 


a 


.(out) 


2V 


V<£ (2 ’ 0) vV 2,0) l - V 2 |V(/? (2 ’ 0) | 


8(3 + 2 u)<j% 

Noting the Laplacian acting on the 2nd term we define 

“ -|V^ 2 ’°)| 2 , 


mix ) = ^ 2 ~ 2S > + 


8(3 + 2w)4° Ut) 
so that our equation after using (3.15) becomes 

8itG 


(3.20) 


(3.21) 


vV mix > = -- 


a 


4(3 + 2c o)^ ut) 3 + 
Then we can solve the above equation to get 


V 


pVip 


( 2 , 0 ) 


(3.22) 


^ mix \t,x) = 


a 


2 G 


d 3 x'-, 


1 


4(3 + 2ui)^° ut) 3 + 2w J I*-*'I 
After some integration by parts and using the identity 

1 1 


V* 


p(t,f')V x ^ ( 2 ’ 0 ) (t, x') . (3.23) 




\x' — x\ 


\x' — x | 3 


(x -x), 


(3.24) 
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we find 




(mix) 


(t,x) = 


G C 

4(3 + 2w)4° ut) V 2 + w 


J o 

X f d?x' ( dV ■ (x' - x"). 

I I \x — x'\ 6 \x — x"r 


(3.25) 


To proceed we also need the term |Vc/^ 2,0 )| 2 which is given by 


a 


8(3 + 2tu)$ ut) 


|V ^ (2 ’ 0) | 2 = 


a 


G c 


8(3+ 2u )](j) { 0 out) V 2 + w 
x [ d 3 x' [ d 3 x" p(x')p(x") 


,3,„/ / ^3_// „/„//>> (g ~ ~ &') 

I I 3 I djQ _ /•£»// I 3 


(3.26) 


so that the full solution is 




(2 ’- 2 >(*,*) = - 


a 


1 2+u I 


(^y J d 3 x' j d 3 x"p(t,x')p(t,x") 


8(3 + 2u)<f> K Q 

(x — x') • (x — x") 2 — •'?') • (a: 7 — x") 


I I ^ I _ cc ^ ^ I ^ 


q™ _ I 3 I /yJ rjj^J t I 3 


(3.27) 


3.2.6 The PPNV metric 

The new potential (3.27) that we have just found, is an example of what we call a Post- 
Newtonian-Vainshteinian potential. To follow the spirit of the PPN formalism, let us define 
it slightly differently, i.e. without the constants that appear in front of the solution. We 
define the PPNV potential as 


Uy Ut ^ = [ d 3 x' [ d 3 x" p(t, x')p(t, x") 


(x — x') ■ (x — x") ^(x — x') ■ ( x 1 — x") 


I _ nj^l I 3 I f I 3 


I 3 I nj^J _ t I 3 


, (3.28) 


so that up-to 0(2, —2) the solution for is 




( 2 ,<|- 2 |) = 2G jj 


a 


f Gc V 

\ 2 +w ) 


-u, 


3 + 2lu~ 8(3 + 2cu)$ ut) V 


(out) 


(3.29) 


where the notation < | — 2n| means to sum all orders less than or equal to | — 2n|, i.e. 

^ 2 ,<|- 2 | = ^ 2,0 + 

We may now use (3.8) and (3.12) to determine the metric up-to 0(2, —2). From (3.8) 
the metric solution to 0(2, —2) is 


h ( 0 2 o’- h2l) = 2 G C U + 2g v G 3 c u\° ut) . 


(3.30) 


The parameter gy is a PPNV parameter that for the Cubic Galileon is 

l 3 


7T 


gv = 


M n 


4 [ (2 + w)Aj ’ 


(3.31) 
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and for Brans-Dicke (GR) gy = 0. Likewise from (3.12) the metric solution to 0(2, —2) is 


hij 


2 7 G C U + 7 V G 3 c ul out) 


7 iji 


(3.32) 


where 7 y is a new PPNV parameter that measures Vainshteinian corrections. In the Cubic 
Galileon theory we have 


7 r 

iv = -gv = ^ 


Mp y 

(2 + w)Aj ’ 


(3.33) 


and in Brans-Dicke or in GR 7 y = 0. Note that as uj — > 00 we recover GR and the Vainshtein 
corrections vanish. This makes sense since in that limit Co — > 00 and the scalar completely de¬ 
couples. Indeed, by canonically normalising we see that the Galileon interaction is suppressed 

3 

by a divergent strong coupling scale A strong ~ Acq . 


3.3 Dualizing the Cubic Galileon 

When the observer is placed close to the source the held self-interactions become strong and 
the standard approach used in the outside region breaks down. In order to proceed we must 
recast the original theory into a classical dual formulation [51, 52]. The basic idea is to 
introduce auxiliary variables “dual” to the interaction terms in the standard action by a 
Legendre transform. It turns out that the equations of motion for the new set of variables 
can be treated perturbatively since the expansion parameter 1/A is flipped to A. 

A technique to compute such a dual action was initially given in [51], but in their dual 
variables the resulting equations involve non-analytic functional forms that make the math¬ 
ematics quite cumbersome. A neater and more direct approach was put forth in [52], where 
instead of using Legendre transforms as dual variables, the interaction terms are directly 
related to the auxiliary variables by using Lagrange multipliers. In these variables simpler 
analytic equations arise keeping all the advantages of the original approach. The dual action 
for the Cubic Galileon is given by 

Sdual = ^/ **Sl{*R- |(V«> +V (^ + »V^) - ■ 

(3.34) 

where A^ is the dual field corresponding to and Z is the dual field for □</>. By extremizing 

the action with respect to the field and after some algebraic manipulation we find the two 
relations between the dual fields and gradients of 4> 


V M </> = a 2 Ap, 

□ (/ = Q~ 2 Z, 


the generalized Einstein equations 


1, 


(frRpv — 8ttG [ T^u n T ) + ct 2 V ^ A u + n Zg^ u + ± ol 2 A^A V 


1 


UJ 


+ 


1 

8^3 


1 


1 A 2 


— ZA /Ji A u + AqjVy )A + —ZA^g^ — 3a 2 ——A^ j A l , 


(3.35) 

(3.36) 


(3.37) 
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and the field equation for the dual fields (which is equivalent to the scalar equation in the 
standard formulation) 


— OA 2 + 2V\(ZA x ) 


— 647rG(/> 2 T + a 2 


8(2w + 3)(j) 2 Z 


b A^^A 2 + + 

9 <P 


_ 1 
9a 2 

<f > 2 



(3.38) 

Once again, we have manipulated the field equations to put them into a form more useful 
for performing the PPN expansion. It can easily be checked that eliminating the dual fields 
from (3.37) and (3.38) using (3.35) and (3.36) recovers the field equations in the standard 
formulation. 

In the dual formulation the expansion parameter is cV 1 / 2 which in this case will have 
order 0(1, 1). 


3.4 The PPNV expansion inside the Vainshtein radius 

As with the outside region, we expand the scalar field as 

<t> = rf n) (! + ¥>) = tf n) (l + <P {2) + V (4) + ■••), (3-39) 

where is a constant and p a perturbation. Note that due to the non-perturbative nature 
of the theory at the Vainshtein surface, the constant will in general not be equal to (j)Q OUt ^. 
However, they can in principle be related by appropriately matching the scalar </> across the 
two regions. 

Using the definitions of the dual fields (3.35) and (3.36) we find that to PPN order 2 

4 in) ViP (2) = a~^A[ 1] = <$ n) ar*ViBW, (3.40) 

rf n) vV 2) = cT5Z (1) = ^ n) Q-sV 2 B (1) , (3.41) 

where we have defined the scalar field B whose gradient gives the dual field A^. Thus, in 
general, the orders of <p will be one greater than the orders of the B field (either PPN or 
Vainshtein order). Notice that as the constant a ~2 is of order 0(1, 1), it increases both the 
PPN order and the Vainshtein order of any terms multiplying it by one. This means that 
the lowest PPN order for the scalar perturbation p is 2 while the lowest Vainshtein order for 
the scalar field perturbation p is 1 (given that the lowest possible Vainshtein order for the B 
field is zero). 

We now consider the generalized Einstein equations to 0(2) in PPN. 


3.4.1 fioo to 0(2) 

Taking the Rqq equation (3.37) and expandng to 0(2) using also (2.5) we find 

Iv 2 4 2 0 } = -4 t TG N p+ l -a-H 2 B^\ 


where we have defined 


G n = 


G 

iO) ‘ 
9 o 


We can determine /iqq completely in terms of the new potential B^. The answer is 

h$ = 2G n U + a~ ^B^. 


(3.42) 

(3.43) 

(3.44) 
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Since the lowest Vainshtein order for B^ is 0 taking the limit a~ 2 —> 0 recovers General 
Relativity. Thus the Vainshtein mechanism is directly manifested in the dual formalism. In 
other words, to 0(2, 0) the B ^ term does not contribute and we recover the usual Newtonian 
potential 

4o 0) = 2 G n U. (3.45) 

Hence Gjsr has the meaning of the observable Newtonian gravitational constant inside the 
Vainshtein radius, i.e. the Newtonian constant measured by a table-top experiment. 

Further out away from the source, the Vainshteinian potential B^ will start to have 
an effect and we need to take into account the Vainshteinian corrections. This is achieved 
via the scalar equation to 0(2) (see Eq. 3.51) 

3.4.2 hij to 0(2) 

Taking the Rij equation (3.37) and expandng to 0(2) using also (2.6) we find 



Vi/T\- + V-ih 


V hij + o VjVj(/lg 0 ^ 


+a 2 


2 J 2 
ViVj-B (1) + Jv 2 R (1) 


7 ij 


h) = AnGNp-fij 


Imposing the gauge fixing 


V k h k i = Vi {^-h- l -h { S+a-^B^, 

turns (3.46) into 

- 2 hij = i'rrGNp'jij + V 2 B^jij, 

so that we obtain hjj as 

hij = (2 G N U-a-^B^ 7ij , 


(3.46) 


(3.47) 


(3.48) 


(3.49) 


We notice that no new potential (apart from B W) arises. Once again as a —> 0 we 

recover GR, i.e. 

h\f 0) = 2G N U lij . (3.50) 

3.4.3 Scalar equation to 0(2) 

Just as in the case outside the Vainshtein radius, in order to proceed further we need the 
scalar equation to 0(2). We find 


V 2 (|VR (1) | 2 ) - 2V • (V 2 R (1) VR (1) ) = -64t xGp - 8cT 1/2 (2u; + 3 )<^ n) V 2 R (1) , (3.51) 

we then expand the dual field in Vainshtein orders, i.e. 

OO 

R(i) = + 5 (1>1) + R {1 ’ 2) + ... = Y^ R (1,n) - (3.52) 

n =0 
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Taking the lowest order we get an equation for the dual field B ^ 1 '°) which is 

V 2 (jVB (1 ’ 0) | 2 - 16 Gu} - 2V • (v 2 R (1 ’ 0) VR (1 ’ 0) ) = 0. (3.53) 

However, at the same limit, the terms involving B ^ disappear in the metric solutions (3.44) 
and (3.49). What happens is that the solution for the metric to Vainshtein order 0 is simply 
General Relativity while the field B W decouples and obeys its own differential equation with 
no bearing on the metric or matter evolution. 

But the field B^ will eventually have a bearing on the metric. Solving (3.53) and 
inserting into (3.44) and (3.49) gives us the first correction to the GR metric, i.e. the 
metric solution up-to 0(2,1). Unfortunately a complete analytic solution to (3.53) is (to 
our knowledge) impossible, except in idealized situations like spherical symmetry. Indeed, 
by simple shuffling of the terms involved, (3.53) can be re-written in the form 

V^B^ViVjB^ - (V 2 R (1 ’ 0) ) 2 = —32vr Gp. (3.54) 

We recognise the above equation as a Monge-Ampere like equation whose general solution 
is unknown, and one has to resort to numerics. Thus, it is impossible to find this new 
Vainshteinian potential as a closed form integral as we have done in the case outside the 
Vainshtein radius. 

3.4.4 Reconstructing the metric: the 0(2) metric solution to all Vainshtein 
orders 

Fortunately, the polynomial structure of (3.51) lends itself to the use of perturbation theory. 
Once the lowest order solution is found, for instance £>1 1,0 ) then higher orders B ( ' 1,n ' > will obey 
a linearized equation coming from (3.51). More specifically, expanding order-by-order (3.51) 
may be re-written as 

£(B^ 1,0 \ r( 1,0 )) = 64vr Gp, 

£(R (1 ’ 1) ,R (1 ’ 0) ) + C(B^'°\B^’ 1) ) = 8a"^(2w + 3)4 in) V 2 R (1 ’ 0) . (3.55) 


The operator £ is the bi-linear operator £ : C°° x C 

C(u,v) = 2V • (V 2 u)Vu - 


00 -> C°° defined by 
- V 2 (Vri- Vv), 


(3.56) 


for two arbitrary functions «,»£ C°°. The operator £ is a non-symmetric (C(u, v) / C(v, u)) 
bi-linear form, i.e. it obeys the properties C(u + v,w) = C(u,w) + C(v,w) and C(\u,v) = 
C(u, \v). Then we can formally write the solutions to the hierarchy of linearizations (3.55) 
by considering the operator £ : C°° C°° defined by 

CC(u,v) = uv. (3.57) 

The operator £ is distributive with respect to addition and to multiplication by a constant 
that ensures that £ is invertible and a solution can be constructed. 
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In the spirit of the PPN potentials and also the Vainshteinian potential Uy Vt> the last 
considaration is to rescale the potentials to factor out any model-dependent parame¬ 

ters. We therefore define the series of Vainshteinian potentials for each Vainshtein order n 
as 


Tj{in) = 




4 y/G |_2(2o; + 3)0Q n) 

so that the hierarchy of linearizations (3.55) becomes 


B^' n \ 


(3.58) 


C(ut\ut ] ) = 4t rp, 

C{u[ in \u^ in) ) + C(U^ n \u[ in) ) = v 2 u^ in \ 


(3.59) 


The above hierarchy can then in principle be solved order-by-order where now the only input 
is the matter density p. Formally this is achieved via the operator £ that ensures that the 
solutions exist: 

u { 0 in) = 

U i n) = ^n)^ V2U 0 in) )^ (3-60) 


Once this is done the metric can then be reconstructed as 

4o = 2G N U + 2Y j g% n) G- N n U^\ 


h-ij = (^G n7 U + 2 7lt n) Gjrut^J 7m 


(3.61) 

(3.62) 


where we have introduced the PPNV parameters (jn"' 1 and 7 along with the PPN param¬ 
eter 7 . For GR we have 7 = 1 and gn^ = 7 = 0, for the Brans-Dicke theory 7 = and 
Qrl^ = 7 = 0 while for the Cubic Galileon we have 7 = 1 and 


ffi m) = 


-tRO = 

in 


2 1+n (2 u + 3) r 


(8 


n +1 

7r) 2 


A V^ 


,(in) 

^0 


3(n+l) 


M (n°m) 


(3.63) 


where Mp lom ' > = is the nomimal Planck mass (defined using the measured Newton’s 

constant Gn rather than G. 


4 Back to spherical symmetry 

We have determined the PPNV expansion for the Cubic Galileon theory to 0(2) in PPN both 
inside and outside the Vainshtein radius in terms of general potentials. Let us now compare 
the expansion with the results for spherical symmetry that we have found in section 2.2. We 
start from the expansion outside the Vainshtein radius. The second term inside the integral of 
the Vainshteinian potential (3.28) is antisymmetric in the exchange of x! and x" which means 
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that it will vanish in any spherically symmetric situation. Using p(t,x) = M8^ 3 \r) = If-yw 
we evaluate the integral as Uy = Since the Newtonian potential is U = we determine 
the metric solution as 


hoo = 


hij — 


2 G C M 
r 

2 G C M 
r 


1 + 


gyGlM 


2 G C M 


i-*->(t) 


1 +oj 'yyG^M 

V *3 


where f3 [out) = 


2 +LU 
1 


7 ij — 


2 G C M 


r 

1 + OJ 


2 + U + P(out) ( r ) 


(4.1) 
7ij> (4-2) 


which is the spherically symmetric solution in section 


64tt(3+2w) 2 (2+w) (^° ut) ) 2 
2.2 (only now we have also found the metric component h tJ . 

We turn now to the inner region. The operator £ for spherical symmetry gives 

2 4 

L{u, v) = u v — u v -\—(uu v . 


(4.3) 


For the Vainshtein order 0 case, u = v = Uq 11 ' 1 so that C{u,u) = ^z-^['r(u') 2 ] = ^5(r) 
which gives the solution £/q = ±\]Mr. Choosing the minus sign of the solution, i.e. 

Uq 71 ' = —yjMr, leads to the metric solution 


>>w = 2 


hij — 2 G N 


M 


1 - 4v / 27r(0Q n) ) 3//4 


1 + 4v/2)((0[ ) in) ) 3 / 4 


(nom) 


(nom) 


3/2' 


3/2' 


7ii> 
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(4.4) 

(4.5) 


(nom) 


M, 


1/3 


where we have defined the nomimal Vainshtein radius ry° m ‘ = 

The plus/minus signs in the solution for Uq"^ has to do with the existence of solution 
branches in the field equations. This is the same issue that we encountered in (2.10) where at 
that instance we chose the plus sign in front of the square root. Had we chosen the opposite 
sign, we would have recovered the solution found here with the plus sign, i.e. Uq"^ = \jMr 
for the inner region r < ry. However, for the outer region, the leading order term would 
have then been x ~ ? ’ 2 > rather than ~ with the constant of proportionality depending only 
on a and not the mass of the source. Therefore this branch of the solution is not present 
in our PPNV expansion in the outer region r > ry for the reason that it is not sourced by 
matter (and therefore cannot be recovered from (3.28)) but also due to our assumption of 
a background Minkowski metric, an assumption which would be violated by this solution 
branch. Perhaps, one has to expand around de Sitter space in order to have this branch in 
the solution space. 

Let us now turn to the values ^ and (j)^ ut \ In principle, one can rescale units, for 
instance, A —»• A\J 0q OU< ^ and similarly for all other mass (and length) scales. Then one can 
set either of <or (but not both) to unity by simple unit re-definition. It is a matter 

of choice which one, but for what follows we choose to set = 1 as it is the most relevant 
to the solar system (in which case Gjy = G and Mp lom ' > = M p ). 
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In analogy with the standard PPN formalism, we may try to check the effect of the 
parameters g n and on null geodesics. Using the geodesic equation for light we find that the 
Shapiro time-delay and the light deflection angle are both proportional to the combinations 
Sn + 7 n for orders n. While these may be used to place constraints on the combination 
9n + In, clearly, the Cubic Galileon terms can not be constrained using light, although a 
disformally coupled Galileon would presumably give a different result. Of course, by tracking 
material bodies, it may also be possible to place constraints on the PPNV parameters. We 
leave these considerations for a future work. 

For completeness, let us give a rough estimate of the deviations from GR in the solar 

( 2 ) 

system. In figure 1 we display the functions h (m and hij versus r as well as the residuals 
from GR. We assume a spherically symmetric solar system with all matter concentrated 
in the sun. Then Mq/M p = 1.12 x 10 5 ' so that ry = 1.02 x 10 -5 au. For a choice of 
A ~ (■ M pl H 2 0 )3 1/1000km ~ 2 x 10 13 eV (ry ~ 5 x 10' au) (in order that the scalar modifies 
the late time cosmology from the present epoch onwards), this theory predicts deviations of 
the order ~ 10 -10 at the orbit of Jupiter and ~ 10 -9 at the orbit of Neptune. It is interesting 
also that the Vainshtein radius of the Sun for this choice of A is larger than the distance to 
the nearest star. 



Figure 1 . Solutions for the potential h t j for the Cubic Galileon inside the Vainshtein radius. We 
display the Newtonian term (dash), the Vainslrteinian correction (dotted), the total potential (solid) 
and h% 1] /h% 0) (dash-dot). 

5 Conclusion 

We have presented a general scheme for calculating the metric for theories exhibiting Vain¬ 
shtein screening, in which there is a unique scale beyond the Schwarzschild radius where 
non-linearities begin to become important. This scheme, the Parametrized Post-Newtonian 
Vainslrteinian expansion, is an expansion in two small parameters, namely, small velocities 
and a parameter that determines the strength of the Vainshteinian corrections. We saw that 
such theories have two regimes: inside the Vainshtein radius and outside. One has to perform 
a separate expansion in each regime. In particular, the inner regime requires the use of a 
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dual and classically equivalent action that inverts the scale A and makes the expansion in a 
small parameter feasible. 

The long term goal of the PPNV framework is to help determine the compatibility of 
the relevant theories with solar system and other strong-field data and we expect it to play a 
fundamental role in constraining them independently from cosmology. This paper should be 
regarded as the first in developing a full formalism. Whilst we have explained the formalism 
for the class of theories under consideration i.e. those with a single Vainshtein radius, we 
have only presented this explicitly for the Cubic Galileon, working up to PPN order 2. We 
expect new PPNV potentials to arise as we introduce more operators and/or increase the 
order of our expansion, just as is the case with PPN. In a separate calculation we will present 
the expansion to 0(4) in PPN for the Cubic Galileon. It would also be pertinent to identify 
PPNV potentials for the higher Galileons, leading order k-essence terms, and for massive 
gravity (to name a few). 

For the Cubic Galileon case, we have explicitly demonstrated the working of the Vain¬ 
shtein mechanism. Interestingly, we found that the theory tends to Brans-Dicke far outside 
the Vainshtein radius, i.e. in the cosmological regime, while it tends to GR deep inside the 
Vainshtein radius, where solar-system tests lie. Thus, the Cubic Galileon is a concrete exam¬ 
ple where solar-system constraints on Brans-Dicke theory do not apply but rather one has 
to use cosmological constraints [56]. 

A more challenging extension would be to allow for multiple macroscopic scales, mul¬ 
tiple Vainshtein radii if you like. This will obviously mean expanding in more than two 
variables. Note that the classical dual can be formulated to cope with mutliple strong cou¬ 
pling scales (which would in turn lead to multiple Vainshtein radii). Difficulties could occur, 
however, when interactions suppressed by the same strong coupling scale kick in at differ¬ 
ent macroscopic scales due to their differing derivative structure. It is not entirely clear 
how the classical dual will capture this as an expansion in the appropriate Vainshtein radii. 
Presumably we will have to extend the classical dual itself in order to cope with this. 

Finally, let us comment on the so called self-screening effect which has been argued to 
play an important role for Vainshtein systems [47, 57], and suggests, for example, that the 
moon cannot be treated as a test body in the field profile of the earth. We essentially agree 
with this. The point is that our method allows us to extract a perturbative solution for the 
scalar profile and the geometry inside and outside the Vainshtein radius for the entire system. 
If we take our probe to be part of the source in that calculation, then self-screening effects 
should be accommodated although to be practical one may wish to include an additional 
expansion is terms of, say, the relative mass scales in a two body problem. The non-linear 
nature of the inner PPNV expansion would be key to incorporating self-screening. In contrast, 
if we do not include our probe as part of the source, then we are genuinely treating it as a 
test particle and ignoring any self-screening. 

To be more explicit, a small body orbiting the Sun can indeed be taken as part of the 
source so that p = Psun+Pbody- We then solve our field equations and calculate the corrections 
to the potential as explained in this paper. Of course there is no spherical symmetry but 
that is a calculational problem rather than an issue with our formalism. Once we have the 
metric, the small body will certainly follow the geodesics of that metric just like a test body. 
The key point, however, is that it will not follow the geodesics of a metric calculated from 
/9 sun alone. This is the correct statement regarding both the self-screening papers and our 
work. How big the effect is, and how one can re-cast it as a violation of the SEP, are two 
different questions, but two answerable questions within our formalism. 
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